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Abstract 

The f(T) theory, which is an extension of teleparallel, or torsion scalar T, gravity, is recently 
proposed to explain the present cosmic accelerating expansion with no need of dark energy. 
In this Letter, we first perform the statefinder analysis and Om(z) diagnostic to two concrete 
f(T) models, i.e., f(T) = a(—T) n and f(T) = —aT(l — e pT °/ T ), and find that a crossing of 
phantom divide line is impossible for both models. This is contrary to an existing result where 
a crossing is claimed for the second model. We, then, study the constraints on them from the 
latest Union 2 Type la Supernova (Sne la) set, the baryonic acoustic oscillation (BAO), and the 
cosmic microwave background (CMB) radiation. Our results show that at the 95% confidence 
level n m0 = 0.272±g;gjjfj, n = 0.04+°; 22 for Model 1 and n m0 = 0.272±g;g||, p = — 0.02±g;|J for 
Model 2. A comparison of these two models with the ACDM by the XMi n /dof (dof: degree 
of freedom) criterion indicates that ACDM is still favored by observations. We also study the 
evolution of the equation of state for the effective dark energy in the theory and find that Sne 
la favors a phantomdike dark energy, while Sne la + BAO + CMB prefers a quintessencedike 
one. 
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I. INTRODUCTION 

Various cosmological observations, including the Type la Supernova [lj], the cosmic 
microwave background radiation Q and the large scale structure [3, 4|, et al., have revealed 
that the universe is undergoing an accelerating expansion and it entered this accelerating 
phase only in the near past. This unexpected observed phenomenon poses one of the 
most puzzling problems in cosmology today. Usually, it is assumed that there exists, in 
our universe, an exotic energy component with negative pressure, named dark energy, 
which dominates the universe and drives it to an accelerating expansion at recent times. 
Many candidates of dark energy have been proposed, such as the cosmological constant, 
quintessence, phantom, quintom as well as the (generalized) Chaplygin gas, and so on. 
However, alternatively, one can take this observed accelerating expansion as a signal of 
the breakdown of our understanding to the laws of gravitation and, thus, a modification 
of the gravity theory is needed. One of the most popular modified gravity models is 
obtained by generalizing the spacetime curvature scalar R in the Einstein-Hilbert action 
in general relativity to a general function of R. The theory so obtained is called as the 
f(R) theory (see |5| for recent review). 

Recently, a new modified gravity by extending the teleparallel gravity |6|] is proposed 
to account for the present accelerating expansion 7H10I]. Differing from general relativity 
using the Levi-Civita connection, in teleparallel gravity, the Weitzenbock connection is 
used. As a result, the spacetime has only torsion and thus is curvature-free. Similar to 
general relativity where the action is a curvature scalar, the action of teleparallel gravity 
is a torsion scalar T. In analogy to the f(R) theory, Bengochea and Ferraro suggested, 
in Ref. p|, a new model, named f(T) theory, by generalizing the action of teleparallel 
gravity, and found that it can explain the observed acceleration of the universe. Let us 
also note here that models based on modified teleparallel gravity may also provide an 
alternative to inflation [l^]. Another advantage the generalized f(T) torsion theory 
has is that its field equations are second order as opposed to the fourth order equations 
of f(R) theory. More recently, Linder proposed two new f(T) models to explain the 
present cosmic accelerating expansion |8| and found that the f(T) theory can unify a 
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number of interesting extensions of gravity beyond general relativity. In this Letter, we 
plan to first perform a statefmder analysis and an Om diagnostic to these models and 
then discuss the constraints on them from the latest observational data, including the 
Type la Supernovae released by the Supernova Cosmology Project Collaboration, the 
baryonic acoustic oscillation from the spectroscopic Sloan Digital Sky Survey, and the 
cosmic microwave background radiation from Wilkinson Microwave Anisotropy Probe 
seven year observation. We find that for both models the crossing of the —1 line is 
impossible. This is consistent with what obtained in Ref. [10|, but in conflict with the 
result obtained in Ref. [8j where a crossing is found for the exponential model. 



II. THE f(T) THEORY 

In this section, following Refs. {7, 8|, we briefly review the f(T) theory. We start with 
teleparallel gravity where the action is the torsion scalar T defined as 

T = S/ V T\ V , (1) 

where T a is the torsion tensor 

r%, = eS(0„e? - d„ei) , (2) 

and 

Sr = \ (K^ a + 5£T™ a - 8»T^ a ) . (3) 

Here is the orthonormal tetrad component, where A is an index running over 0, 1, 2, 3 
for the tangent space of the manifold, while /i, also running over 0, 1, 2, 3, is the coordinate 
index on the manifold. The spacetime metric is related to through 

= riABefe^ , (4) 

and K^ p a is the contorsion tensor given by 

K^ a = - l -(T^-T^-T a n (5) 
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By assuming a flat homogeneous and isotropic Friedmann-Robertson- Walker universe 
which is described by the metric 

ds 2 = dt 2 - a 2 {t)5 ij dx i dx j , (6) 

where a is the scale factor, one has, from Eq. (CQ), 

T = -6H 2 , (7) 

with H = da -1 being the Hubble parameter. 

In order to explain the late time cosmic accelerating expansion without the need of 
dark energy, Linder, following Ref. [/J, generalized the Lagrangian density in teleparallel 
gravity by promoting T to be T + f(T). The modified Friedmann equation then becomes 

H 2 ~P~{- *H 2 f T , (8) 

2 16nGP + 6H 2 + f + 12H 2 f T 
1 ' 2 AH 2 f xt - 2 - 2f T ' 1 ' 

where a prime denotes a derivative with respect to In a, p is energy density and P is the 

pressure. Here we assume that the energy component in the universe is only matter with 

radiation neglected, thus P = 0. 

From Eqs. (JHHJJ), we can define an effective dark energy, whose energy density and the 

equation of state can be expressed, respectively, as 

^"-itd ( - ! + 2 ™ (10) 

f/T - fr + 2Tf TT 
Weff ~ "(l + /r + 2T/ Tr )(//T-2/ r ) ' (U) 

Some models are proposed in Refs. [l, §| to explain the present cosmic accelerating 
expansion, which satisfy the usual condition f/T — > at the high redshift in order 
to be consistent with the primordial nucleosynthesis and cosmic microwave background 
constraints. Here we consider two models proposed by Linder 

• Model 1 



f(T) = a(-T) 
4 



(12) 



Here a and n are two model parameters. Using the modified Friedmann equation, one 
can obtain 

01 = (6 ^ )1_ni 2^f ' (13) 
where f2 m o = 87r ^°- > is the dimensionless matter density today. Substituting above ex- 
pression into the modified Friedmann equation and defining E 2 = H 2 /Hq, one has 

e 2 (z) = a m0 (i + z f + (i - n rn0 )E 2n . (14) 

Let us note that this model has the same background evolution equation as some phe- 
nomenological models [3, [l4| and it reduces to the ACDM model when n = 0, and to 



the DGP model 



15| when n — 1/2. When n — 1, the Friedmann equation (Eq. (jSJ)) can 



be rewritten as H 2 = „f? G , p, which is the same as that of a standard cold dark matter 

3(1— a) "> 

(SCDM) model if we rescale the Newton's constant as G — > G/(l — a). Therefore, in 
order to obtain an accelerating expansion, it is required that n < 1. 
• Model 2 

f(T) = -aT(l - e pTo/T ) , (15) 

which is similar to a f{R) model where an exponential dependence on the curvature scalar 
is proposed [if], Q • Using the modified Friedmann equation again, we have 

a= l-(l-2p)eP' (16) 

and 

E 2 - E 2 e p / E2 4- 2ne p / E2 

e 2 {z) = n m0 (i + zf + (i - n rn0 r 1 _( 1 _ 2 ^ )eP • ( 17 ) 

It is easy to see that p = corresponds to the case of ACDM. 



III. STATEFINDER ANALYSIS AND Om DIAGNOSTIC 

In order to discriminate different dark energy models from each other, Sanhi et al. 
proposed a geometrical diagnostic method by adding higher derivatives of the scale fac- 
tor 



181 ] . In this method, two parameters (r, s), named statefinder parameters, are used, 



which are defined, respectively, as 



r = 



aH 3 ' 

r - 1 



S = 



(19) 



3(g-l/2) ' 

where q = — -jp- is the decelerating parameter. Apparently, ACDM model corresponds to 
a point (1, 0) in (r, s) phase space. The statefinder diagnostic can discriminate different 



models. For example, it can distinguish quintom from other dark energy mode 

The Om(z) is a new diagnostic of dark energy proposed by Sahni et al. 20]. It is 
defined as 



19]. 



Om(z) = 



EHz) 



1 



(20) 



Apparently, this diagnostic only depends on the first derivative of the luminosity D^(z) 
(see Eq. (1A3j) ). Thus, its advantage, as opposed to the equation of state of dark energy, is 
that it is less sensitive to the observational errors and the present matter energy density 
Qmo- One can use this diagnostic to discriminate different dark energy models by exam- 
ining the slope of Om(z) even if the value of Q m o is not exactly known, since the positive, 
null, or negative slopes correspond to w < —l,w = — 1 or w > —1, respectively. 
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FIG. 1: The evolutionary curves of statefinder pair (r, s) (left), pair (r, q) (middle) and Om{z) 
(right) for Model 1 with Q m0 = 0.278. 



Here, we perform the statefinder and Om diagnostics to two f(T) models, i.e., Model 
1 and Model 2 given in the previous section. In Figs. ([1]) and ([2]), we show the diagnostic 
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FIG. 2: The evolutionary curves of statefinder pair (r, s) (left), pair (r, q) (middle) and Om{z) 
(right) for Model 2 with fl m0 = 0.278. 

results with Q m0 = 0.278, which is the best fit value obtained from Sne la and BAO with 
a model independent method 2l| and is also consistent with the result in the next section 
of the present Letter. The left panels show the evolutionary curves of statefinder pair 
(r, s), the middle panels are the evolutionary curves of pair (r, q), and the right panels are 
the Om(z) diagnostic. Although, both Model 1 and Model 2 evolve from the SCDM to a 
de Sitter (dS) phase as one can see from the middle panels of these figures, the effective 
dark energy for Model 2 with p ^ is similar to a cosmological constant both in the high 
redshift regimes and in the future, while for Model 1 with n ^ this similarity occurs 
only in the future. 

As demonstrated in Ref. 22], for a simple power law evolution of the scale factor 
a{t) ~ t 2 / 3 " 7 , one has r = (1 — 37)(1 — 37/2) and 5 = 7. Thus, a phantom-like dark 
energy corresponds to s < 0, a quintessence-like dark energy to s > 0, and an evolution 
from phantom to quintessence or inverse is given by a crossing of the point (1,0) in (r, s) 
phase plane. A crossing of phantom divide line is also represented by a crossing of the 
red solid line (ACDM) in middle panels ((r, q) plane) of Figs. ([HE]). Therefore, we find, 
from the left and middle panels of Figs. (pQ 12]), that n > (Model 1) or p < (Model 
2) f(T) corresponds to a quintessence-like dark energy model, while n < (Model 1) or 
p > (Model 2) corresponds to a phantom-like one. A crossing of the phantom divide 
line is impossible for Model 1 and Model 2. These results are also confirmed by the 



Om(z) analysis given in the right panels. In order to further confirm our results, we 
redo our analysis with other values of fi m o, such as Q m o = 0.2 or 0.5, and find that the 
result remains unchanged. Thus, we conclude that the phantom divide line is not crossed 
for both models. This is in conflict with what given in Ref. 8] where a crossing of the 
phantom line is found for Model 2. 

IV. OBSERVATIONAL CONSTRAINTS 

The constraints on model parameters of Model 1 and Model 2 will be discussed, re- 
spectively, in this section. Three different kinds of observational data, i.e., the Type la 
supernovae (Sne la), the baryonic acoustic oscillation (BAO) from the spectroscopic Sloan 
Digital Sky Survey (SDSS) and the cosmic microwave background (CMB) radiation from 
Wilkinson Microwave Anisotropy Probe (WMAP), will be used in order to break the 
degeneracy between the model parameters. The fitting methods are summarized in the 
Appendix. 

For the Sne la data, we use the Union 2 compilation released by the Supernova Cos- 
mology Project collaboration recently 23]. Calculating the Xs ne i we find that, for Model 
1, the best fit values occur at Q m o = 0.302, n = —0.18 with x\ Hn = 543.953, whereas, for 
Model 2, Q m0 = 0.279, p = 0.08 with Xm n = 543.369. 

Then, we consider the constraints from the BAO data. The_parameter A given by the 
BAO peak in the distribution of SDSS luminous red galaxies [4| is used. The constraints 
from Sne Ia+BAO are given by minimizing Xsne+XBAo- The results are n m0 = 0.279±q;q47, 
n = -0.0lt°,;g (at the 95% confidence level) with xlun = 542.978 for Model 1 and 
tt m0 = 0.278+°;°™, p = 0.02+°;^ (at the 95% confidence level) with xlun = 543.383 for 
Model 2. The contour diagrams are shown in Fig. ([3]). 



R 



urthermore, the CMB data is added in our analysis. The CMB shift parameter 



26 



(all 



27j | is used. The constraints from Sne la + BAO + CMB are given by \\ 
Xsne Jr XBAO Jr XcMB- Fig. (HJ) shows the results. We find that, at the 95% confidence level, 
tt m0 = 0.272+™!, n = 0.04±gJ§ with xliin = 543.168 for Model 1 and VL m0 = 0.272lg;g^, 
p = -0.02i{j:M with xmn = 543.631 for Model 2. 
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FIG. 3: The constraints on Model 1 (left) and Model 2 (right) from Sne la + BAO. The red 
and blue+red regions correspond to 1 — a and 2 — a confidence regions, respectively. 
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FIG. 4: The constraints on Model 1 (left) and Model 2 (right) from Sne la + BAO + CMB. 
The red and blue+red regions correspond to 1 — a and 2 — a confidence regions, respectively. 



With the observational data considered above, we also discuss the constraints on the 
ACDM and the results are fi m0 = 0.277±g;g^| with xlnn = 543.400 (Sne la + BAO) and 
tt m0 = 0.276+^| with xltin = 543.745 (Sne la + BAO + CMB) at the 95% confidence 
level. A summary of constraint results on Model 1, Model 2 and ACDM is given in Table 
(1). From Figs. ()3j HJ) and Table (1), one can see that the ACDM (corresponding to 
n = for Model 1 and p = for Model 2) is consistence with the observations at the 
68% confidence level, while the DGP model (corresponds to n = 1/2 for Model 1) is ruled 
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out at the 95% confidence level. Meanwhile, using the XMin/dof (dof: degree of freedom) 
criterion, we find that the ACDM is favored by observations. 

TABLE I: Summary of the constraint on model parameters and XMin/d°f ■ In the table S+B+C 
represents Sne la + BAO + CMB. 





Model 1 


Model 2 


ACDM 




n 


XMiJ d of 




P 


xlnJdof 




X M iJ d °f 


Sne+BAO 


n 97q +0.050 
u - z ' y -0.047 


-n m +0 - 31 

U.U1_ 54 


0.974 


fl 97O+0.050 
U - 2 '°-0.045 


02 +0 - 48 
U.UZ_ 24 


0.975 


„ 977 +0.040 
u - z ' '-0.038 


0.973 


S+B+C 


n 979+0.036 
u - z ' z -0.032 


n f)4 +0 - 22 

u - u ^-0.33 


0.975 


n 979+0.036 
u - z ' z -0.034 


-0 02+ - 31 


0.976 


2 76+ - 032 

U.ZIU_q 03 g 


0.974 



In addition, we study the evolution of the equation of state for the effective dark 
energy. The results are shown in Fig. (JSJ). The dashed, dotdashed and solid lines show 
the evolutionary curves with the model parameters at the best fit values from Sne la, Sne 
Ia+BAO, and Sne Ia+BAO+CMB, respectively. Apparently, Sne la favors a phantomdike 
dark energy, while Sne la + BAO + CMB favor a quintessence-like one. 
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FIG. 5: The evolutionary curves of the equation of state for the effective dark energy from 
Modell (left) and Model2 (right). The model parameters are set at the best fit values. The 
dashed, dotdashed and solid lines correspond to the constraints from Sne la, Sne la + BAO, 
and Sne la + BAO + CMB, respectively. 
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V. CONCLUSION 



Recently, the f(T) gravity theory is proposed to explain the present cosmic acceler- 
ating expansion without the need of dark energy. In this Letter, we discuss firstly the 
statefinder geometrical analysis and Om(z) diagnostic to the /(T) gravity. Two concrete 
f(T) models proposed by Linder 8] are studied. From the Om(z) diagnostic and the 
phase space analysis of the statefinder parameters (r, s) and pair (s,p), we find that, for 
both Model 1 and Model 2, a crossing of the phantom divide line is impossible, which 



conflicts with the result obtained in Ref. 



where a crossing is found for Model 2. We 



then consider the constraints on Model 1 and Model 2 from the latest Union 2 Type la 
Supernova set released by the Supernova Cosmology Project collaboration, the baryonic 
acoustic oscillation observation from the spectroscopic Sloan Digital Sky Survey Data 
Release galaxy sample, and the cosmic microwave background radiation observation from 
the seven-year Wilkinson Microwave Anisotropy Probe result. We find that at the 95% 
confidence level, for Model 1, Vt m0 = 0.272+°;°!, n = 0.04+°;! with xln n = 543.168 and 
for Model 2, fi m0 = 0.272+°;°^, p = -0.02i°;^ with xlnn = 543.631. We also find that 
the ACDM (corresponds to n = for Model 1 and p = for Model 2) is consistence 
with observations at 1 — o confidence level and it is favored by observation through the 
XMin/dof (dof: degree of freedom) criterion. However, the DGP model, which corre- 
sponds to n — 1/2 for Model 1, is ruled out by observations at the 95% confidence level. 
Finally, we study the evolution of the equation of state for the effective dark energy in 
the f{T) theory. Our results show that Sne la favors a phantom-like dark energy, while 
Sne la + BAO + CMB prefers a quintessence-like one. The analysis of the current pa- 
per also indicates that the f(T) theory can give the same background evolution as other 
models such as ACDM, although they have completely different theoretical basis. Thus, 
it remains interesting to study other aspects of f(T) theory, such as the matter density 
growth, which may help us distinguish it from other gravity theories. 
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Appendix A: Data and Fitting Method 



1. Type la Supernovae 



Recently, the Supernova Cosmology Project collaboration 23J] released the Union2 
compilation, which consists of 557 Sne la data points and is the largest published and 
spectroscopically confirmed Sna la sample today. We use it to constrain the theoretical 
models in this paper. The results can be obtained by minimizing the x 2 value of the 
distance moduli 

XSne - Z2 ' ^ ) 

i=l u,i 

where a 2 i are the errors due to the flux uncertainties, intrinsic dispersion of Sne la absolute 
magnitude and peculiar velocity dispersion. fi t, s is the observed distance moduli and \i t h 
is the theoretical one, which is defined as 

IHh = 5 log 10 D L - fi . (A2) 

Here /i = 5log 10 h + 42.38, h = H / 100 km/s/Mpc, and the luminosity distance D L can 
be calculated by 

D ^ il + ^W)- (A3) 

with E(z) given in Eqs. fTl4| fT7|) . Since fio (or h) is a nuisance parameter, we marginalize 



over it by an effective approach given in Ref. [24j. Expanding Xsne t° Xsne(^o) — An 
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2Bhq + C with A = ]T \ja 2 ui , B = ^[/W 2 ^) - 5 lo Sio Aj/oi* and c = Et/^Oi) - 
5 log 10 D L ] 2 /a 2 j, one can find that x| ne nas a minimum value at /x = B/A, which is given 
by 

B 2 

X.Sne G a ■ 

Thus, we can minimize Xsne instead of x| ne ^° obtain constraints from Sne la. 



(A4) 



2. Baryon Acoustic Oscillation 

For BAO data, the parameter A given by the BAO peak in the distribution of SDSS 
luminous red galaxies |4| is used. The results can be obtained by calculating: 



Xbao 



[A -A 



obs\ 



where A b s = 0.469(n s /0.98) 35 ± 0.017 with the scalar spectral index n s 
the WMAP 7-year data [25[] and the theoretical value A is defined as 



A = n l J n lE{ Zb ) 



4-1/3 



dz' 



z b Jo E(z') 



2/3 



(A5) 
0.963 from 

(A6) 



with zj, = 0.35. 



3. Cosmic Microwave Background 



Since the CMB shift parameter R [26|, |27J contains the main information of the ob- 
servations of the CMB, it is used in our analysis. The WMAP7 data gives the observed 



value of R to be R obs = 1.725 ± 0.018 

as 

R = n 



251 ] . The corresponding theoretical value is defined 



1/2 
mO 



ZCMB dz > 



E Z 



(A7) 



where zcmb = 1091.3. Therefore, the constraints on model parameters can be obtained by 
fitting the observed value with the corresponding theoretical one of parameter R through 
the following expression 



Xcmb 



[R-R, 



obs\ 



a 



(A8) 



R 
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